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Abstract. Olver and Rosenau studied group-invariant solutions of (generally nonlinear) 
partial differential equations through the imposition of a side condition. We apply a similar 
idea to the special case of finite-dimensional Hamiltonian systems, namely Hamilton-Jacobi, 
Helmholtz and time-independent Schrodinger equations with potential on ./V-dimensional 
Riemannian and pseudo-Riemannian manifolds, but with a linear side condition, where 
more structure is available. We show that the requirement of N— 1 commuting second-order 
symmetry operators, modulo a second-order linear side condition corresponds to nonregular 
separation of variables in an orthogonal coordinate system, characterized by a generalized 
Stackel matrix. The coordinates and solutions obtainable through true nonregular separation 
are distinct from those arising through regular separation of variables. We develop the theory 
for these systems and provide examples. 
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1 Introduction 

The primary motivation for this paper was the construction by Olver and Rosenau of group- 
invariant solutions of (generally nonlinear) partial differential equations through the imposition 
of a side condition [27, 28]. Various forms of conditional symmetry methods, including nonclas- 
sical symmetry reduction, have been applied by several authors, e.g. [1, 5, 6, 10, 21, 26, 29]. We 
are interested in using linear second-order side conditions and separation of variables to find 
explicit solutions of Hamilton-Jacobi, Helmholtz, Laplace, wave and heat equations. We exploit 
the special properties of finite order classical Hamiltonian systems and their quantum analogues 
represented by Schrodinger equations to obtain new results on separation of variables. The 
present paper is devoted to Hamilton-Jacobi and Helmholtz or time independent Schrodinger 
equations with potential but the ideas are clearly applicable to more general Hamiltonian sys- 
tems and to non-Hamiltonian systems such as diffusion equations that share some features with 
time-dependent Schrodinger equations. 
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An example is the Schrodinger equation H^f = E^f where H = A + V(x) and A is the 
Laplace-Beltrami operator on some Riemannian or pseudo-Riemannian manifold. We look for 
solutions of this equation that also satisfy a side condition = where S is some given linear 
partial differential operator in the variables x. A consistency condition for the existence of 
nontrivial solutions is [H, S] = AS for some linear partial differential operator A. Moreover 
a linear differential operator L will be a symmetry operator for H, (modulo S*$> = 0) if [H, L] = 
BS for some linear partial differential operator B. To make contact with solutions that are 
separable in some system of coordinates, we restrict to the case where the symmetries and the 
side condition are second-order partial differential operators. 

A Hamilton-Jacobi analog is the equation Ti = E where, in orthogonal coordinates, 

N 

H = J29 jj (x)Pj+V(x.), 

3=1 

which is double the classical Hamiltonian. We look for solutions u(x) of the Hamilton-Jacobi 
equation, where pj = d x jU, subject to the side condition <S(x, p) = 0. The consistency require- 
ment is the Poisson bracket relation {Ti, S} = AS for some function A on phase space. Moreover 
a phase space function C will be a constant of the motion for Ti, (modulo S = 0) if {Ti, £} = BS 
for some phase space function B. 

Finding such systems directly from their definition leads to great computational complexity. 
We explore a new method, based on a generalization of Stackel form, i.e., a generalization of 
separable systems corresponding tho a Stackel matrix, that allows us to generate such systems 
efficiently. 

The second motivation for this paper is the general theory of separation of variables for 
both linear and nonlinear partial differential equations [17, 20, 23, 24, 25]. In these works 
the authors point out that there are two types of variable separation: regular and nonregular. 
Regular separation is the most familiar and was exploited by pioneers such as Stackel [33] 
and Eisenhart [12, 14], [13, Appendix 13]. For regular orthogonal separation of a Helmholtz or 
Schrodinger equation on an A^-dimensional manifold there are always N separation constants and 
the associated separable solutions form a basis for the solution space. There is a well developed 
theory for regular orthogonal separation of these equations, including classification of possible 
separable coordinate systems in various constant curvature spaces and intrinsic characterizations 
of the separable systems, see for example [3, 15, 16, 18, 19, 22, 31, 32] in addition to earlier cited 
references. 

For nonregular separation, on the other hand, the number of separation constants is strictly 
less than N and the separable solutions do not form a basis. Symmetry adapted solutions of 
partial differential equations are prominent examples of this class, but except for these special 
solutions there is virtually no structure or classification theory. First attempts of geometric 
interpretation of nonregular separation are given in [2, 8, 9, 11] where nonregular separation 
is considered as separation in which the separated solution must satisfy additional constraints 
that can be seen as side conditions for the equation. In this paper, however, we will show that 
solutions of Helmholtz and Schrodinger equations with second-order side conditions provide 
a class of nonregular orthogonal separation of variables that can be characterized intrinsically. 
Further, each of these systems is associated with a generalized Stackel matrix and this association 
enables us to generate nonregular separable systems very easily. 

In Sections 2 and 3 we review the Stackel construction for regular additive separation of 
Hamilton-Jacobi equations and regular multiplicative separation and -R-separation for Helmholtz 
equations, and make some comments on their geometric characterizations. We also discuss the 
effect of adding vector and scalar potentials. In Section 4 we introduce a generalized Stackel 
matrix with one arbitrary column and show that its use leads to additive separable solutions 
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of the Hamilton-Jacobi equation, with a side condition. We express the results in a Hamil- 
tonian formalism. Then in Section 5 we carry out the analogous construction for Helmholtz 
and Schrodinger eigenvalue equations. Section 7 is the main theoretical contribution of our 
paper. We show that the requirement of maximal nonregular separation for Hamilton-Jacobi 
equations is equivalent to separation with a generalized Stackel matrix and can be characterized 
geometrically. With some modifications, the same is true of i?-separation for Helmholtz and 
Schrodinger equations (though there is still a "generalized Rodrigues form" gap in the geometri- 
cal characterization). Section 8 is devoted to examples of nonregular separation and discussion 
of their various types and significance. We prove a "no go" theorem to the effect that nonregular 
i?-separation does not occur for the Helmholtz (or Schrodinger) equation with no potential or 
scalar potential on a 2D Riemannian or pseudo-Riemannian manifold. However 2D nonregular 
separation can occur for equations with vector or magnetic potentials. Section 8.3 develops the 
theory for two-dimensional systems with vector potential. The self adjoint Schrodinger equation 
for a charged particle in two spatial dimensions, interacting with a classical electromagnetic 
field, again has no new separable coordinate systems that are obtainable from a generalised 
Stackel matrix. However, nonself-adjoint equations such as the analogous solute transport equa- 
tion with first-order convective terms replacing magnetic potential terms, do indeed have new 
nonregular separable coordinate systems. 

We provide examples of nonregular i?-separation for various zero potential, scalar potential 
and vector potential 3D systems: Euclidean space, Minkowski space and nonzero constant cur- 
vature space, including a Euclidean space example due to Sym [34] of nonregular separation with 
two side conditions. The final Section 9 sums up our conclusions and points the way for future 
research on nonregular separation. The coordinate systems and solutions for true nonregular 
separation, i.e., nonregular separation for which regular separation doesn't occur, are distinct 
from those for regular separation and further study of their scope and significance is in order. 

2 Review of regular orthogonal separation 
for the Hamilton-Jacobi equation 1-i = E 

Write the Hamilton-Jacobi equation in terms of it, = diu(x), as 

N 

n = J2 H r 2 ^ + v^) = E. (2.i) 

1=1 

Here, the metric in the orthogonal coordinates x % is ds 2 = Yl H K dxi f- We want to obtain 

i=l 

additive separation, so that djUi = djdiU = for i ^ j. Requiring that the solution u depends 
on n parameters (Ai, . . . , \n) implies the existence of separation equations in the form 

N 

u 2 + v i (x i ) +Y,s ij {x i )X j = 0, i = l,...,N, \x = -E. (2.2) 

i=i 

Here dkSij(x l ) = for k ^ i and det(sjj) ^ 0. We say that S = (sij) is a Stackel matrix. Set 
T = S~ 1 . 

Then (2.1) can be recovered from (2.2) provided HJ 2 = T lj and V = YjVjT 1 ^. The 
quadratic forms C e = Yl T ^(u 2 + Vj) satisfy C e = —Xg for a separable solution. Furthermore, 

3=1 

setting m = Pi, we claim 

{£?,£>} =0, t + j (2.3) 
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N 

where {H, K,} = {d x iK,d Pi % — d x x%d Pi JC) is the Poisson Bracket. Thus the C , 2 < i < N, are 
i=i 

constants of the motion for the Hamiltonian % = For the proof of (2.3) one notes that 



N 



Differentiating this identity with respect to x l , we find 

N 



J2d i T^s jk {x^+T e % k {x i )=0, 

3=1 



so 

N 



Q.rpij = _ T HJ2s' lk T kj . 



k=l 

We substitute this expression into the left hand side of (2.3) and obtain the desired result after 
a routine computation. 

3 Review of the Stackel procedure for the Helmholtz 
or Schrodinger equation 

We can perform an analogous construction of eigenfunctions for a Helmholtz operator, using the 

N 

Stackel matrix S. We demand eigenfunctions of H in the separated form ^ = f\ ^^{x 3 ') and 

3=1 

depending on the maximal number of parameters. Then, the separation equations are of the 
form 

+ h(x e )d e * + ^(x 1 ) - J28 tj (a*)X^j = 0, £ = 1,...,N (3.1) 
for suitable functions fa, vg to be determined. Thus we have the eigenvalue equations 

N 



L k V = T ki {dj + fade + v e )V = \ k V, k = 1, . . . , N, 



where L\ = H. Based on our calculations of the preceding section, we can establish the com- 
mutation relations 

[L s ,L t ]=0. 

More generally we can consider i?-separation for a general Helmholtz equation. In local coor- 
dinates z 3 on an iV-dimensional pseudo-Riemannian manifold this equation takes the invariant 
form 

H@= | A N + J2F j dj + V ] @ = E&, (3.2) 
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where 

1 N 

V9 i,k=i 

is the Laplace-Beltrami operator. We say that this equation is R-separable in local orthogonal 
coordinates x 3 if there is a fixed nonzero function -R(x) such that (3.2) admits solutions 

N 

G = exp(i?)tf = exp(i?) ]J ¥ j) (x j ), 

i=i 

where ^ is a regular separated solution, i.e., it satisfies the separation equations (3.1). In this 
case the symmetry operators are = exp(i?)L^ exp(— R) and equations (3.1) become 

df 6 + {h{x £ ) - 2d £ R)d t Q + ^v e (x e ) - d u R + (8 e R) 2 - %M A ij 9 = °> 

£ = 1, . . . , N. Then we have 
[L s ,L t } = 0, 

where L\ = H. 

Now consider the case 

HB = {A N + V)Q = EG, (3.3) 

i.e., the case where there is no magnetic field, and V is real. We can define an inner product on 
the space of C°° real valued functions /W(z), f( 2 \z) with compact support in M N , with respect 
to which H is formally self-adjoint 

(/ (1) ,/ (2) )= / f^(z)f^(z)^{z)dz. 



N . N 

If L = Yl ajfc ( z ) , 9?fc + Yl he(z)d£ + W(z) is a real symmetry operator then it can be uniquely 

j,k=i i=\ 

decomposed as L = LW + where is formally self-adjoint and is formally skew- 
adjoint 

1 N 1 N / 1 \ 

v 9 j,k=l V 9 1=1 ^ ' 

Moreover, both L^ l \ are symmetry operators. Note that we have = unless H admits 
a first-order symmetry operator. 

Suppose the system admits N algebraically independent commuting symmetry operators L s 
such that the coefficients (ai k ) of the second-order terms in the symmetries admit a basis 
of common eigenforms. (Without loss of generality we can restrict to the self-adjoint case 
L s = Ls.) Then it is well established [19], from the examination of the third-order terms in 
the relations [L s ,Lt] = that there is an orthogonal coordinate system x J and corresponding 
Stackel matrix [sjk(x 3 )) such that 

N 1 

H = ^dj {hHfdj) + V, h= H X H 2 ■H N . 
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Since the metric is in Stackel form, it is straightforward to verify that the symmetries can be 
rewritten as 




making it evident that the first derivative terms are a gradient. (Here, S is the determinant of 
the Stackel matrix.) Thus via an i?-transform we can express our system in the form 

N 

Then the third derivative terms in the commutation relations [L s ,Lt] = are unchanged and 
the cancellation of second derivative terms is satisfied identically. The first derivative terms just 
tell us that the transformed potential V is a Stackel multiplier, so that it permits separation 
in the coordinates x 3 . The zero-th order relation is satisfied identically. Thus the integrable 
system (3.3) is inseparable. Under the same assumptions, but with a magnetic term added, this 
is no longer necessarily true. It is easy to see that if there is a function G such that F ] = djG, 
i.e., if the magnetic potential is a gradient, then the system is again inseparable. However, if 
the magnetic potential is not a gradient then it is no longer necessarily true that integrability 
implies inseparability, even though the second-order terms in the Laplacian admit a common 
basis of eigenforms. See [4] for some examples. 



4 A generalization of Stackel form 

We define a N x N generalized Stackel matrix by 

/ sn(s 1 ) s 12 (x 1 ) ••• si^ N - 1 (x 1 ) ai(x)\ 
s _ S2i(^ 2 ) s 2 2{x 2 ) ■■■ s 2 ,n-i{x 2 ) a 2 (x) ,^ _ 

\s N1 (x N ) s N2 (x N ) ■■■ s N)N -i(x N ) a N (x)J 

where the Gtj are arbitrary analytic functions of the variables x 1 , . . . ,x . We require that S is 
a nonsingular matrix. Set T = S . Now we assume existence of separation equations in the 
form 

N-l 

u 2 i +v i (x i ) + J2 s id xi )k = °> i = l,...,N, Xi = -E. (4.2) 

(Here Latin indices take values 1, . . . , N and Greek indices take values 1, . . . , N — 1.) Note that 
the term with A at is missing. Thus the general separated solution u will depend on N parameters 
(rather than N + 1), an example of (maximal) nonregular separation. Note that equations (4.2) 
can be considered as the restriction to the case Xjy = of 

N-l 

uf + Vi{x l ) + ^2 s^(x l )X^ + Aatg^x) = 0, 
5=1 

which are not separated for Atv = 0. However, any solution of them is a solution of the N 
equations 

N 

C' = Y, T l i [u] + Vj ) =X e , 1<£<N, C N = 0, (4.3) 

3=1 
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where % = C l is the Hamiltonian. Hence, by solving (4.2) we get a separated solution of 
% = E satisfying also C N = as a side condition. This construction shows that separation with 
a generalized Stackel matrix means (nonregular) separation with a side condition. In this case, 
functions (4.3) become our "restricted constants of the motion" C a = X a (modulo C N = 0). 
Since 

N 

3=1 

differentiating this identity with respect to x l , gives 

N 

2>T« sjtM + T* 4(**)=0, 

3=1 

and 

JV 

Y,(^T £j 8 jN + T t *d i a j )=0, 
3=1 



so 



N-l N 

8iT £j + T ei + T Nj ^ T e %a h = 0. (4.4) 

5=1 h=l 

Using this result it is a straightforward computation to verify the Poisson bracket relations: 

C N . (4.5) 



\k,h=l J 



Relations (4.5) can be considered as the consistency conditions that guarantee the C l are con- 
stants of the motion for the Hamiltonian H, modulo the side condition C N = 0. We have verified 
the relations 

{£\^'} £JV=0 = 0, i,j = l,...,N 

for C 1 = Ti and linearly independent quadratic forms £\ Thus our construction has shown that 
separation with a generalized Stackel matrix implies the existence of N independent constants 
of motion in involution (modulo the side condition) diagonalised in the separable coordinates. 

We can generalize Eisenhart's treatment of Stackel form in which he represented the quadratic 
forms T h in terms of their eigenvalues with respect to the metric T lj = HJ 2 : = p^pHJ^. 

Here, pf ] = 1. Then (4.4) can be rewritten as a system of partial differential equations for 

the pf: 

q tt-2 N-l N 

OiPj + 2 + Pi H i 1^ S itPj + Pj Z> Ph H h d i a h ~ U - ( 4 .b) 

i 5=1 h=l 

In the special case I = 1 these equations reduce to 

o rr-2 AT-1 JV 

, rr-2 / It) , W) rr_ 2 q _ _ n 

^-2 + H i 1^ S i(;Pj + Pj 1^ H h d i a h ~ U - 



3 5=1 /i=l 
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Substituting this result in (4.6) we obtain 



diH: 



N 



*f + ("f - of) i^r + S ("!." - of) >n*°0» = o. 



(4.7) 



?t=l 



Note that if c^a^ = for i ^ h then we recover Stackel form and (4.7) simplifies to Eisenhart's 
equation [12, 14] 



diPf + (pf ~ pT) 



0. 



A way of expressing the identity (4.7) that does not require the introduction of the terms ah 
is to note that at least one of the pi , must be nonzero, say for j = 1. Setting j = 1 in (4.7) 
we obtain 



N 



i 



AN) 



d iP ? + W - p\ 



{(.) {£)\ diH i 



h=i Pi 
Substituting this result back into (4.7) we conclude that 



h; 



d iP ) + {p) - p\ ; )"^ = -jm ( Q iPi + [Pi ~ Pi )-^. 



2 - (N) 
Pi 



-2 



(4.8) 



For future use we remark that if we restrict to the case £ = N then, for Bj = pj/p^, 
(4.8) becomes 



diBj = (Bi - Bj) 9 ^ + Bj(l - Bi) 9 ^, i,j = l,...,N. 



(4.9) 



Remark 1. There is an equivalence relation obeyed by generalized Stackel matrices. If S is the 
matrix (4.1), then for any nonzero function /(x), the generalized Stackel matrix 



/ s 11 (x 1 ) suU 1 ) 



S' 



s^at-i^ 1 ) ai(x)/(x)\ i 
S2,n-i(x 2 ) a 2 (x)/(x; 



S2l(x 2 ) S 2 2{x 2 ) 

\SNi(x N ) s N2 (x N ) ■■■ s N)N -i(x N ) ajv(x)/(x)/ 



(4.10) 



defines exactly the same Hamilton-Jacobi equation, separation equations and side condition as 
does S. 

Remark 2. This construction of Hamilton-Jacobi systems with a side condition can easily be 
extended to construct systems with two or more side conditions. For example, with two side 
conditions C 1 ^^ 1 = 0, C N = 0, the last two columns of the N x N generalized Stackel matrix 
would be arbitrary and the symmetries would be modulo the side conditions: 



*N-1 



i,j = l,...,N. 



In a similar fashion nonregular separability of Helmholtz equations with multiple linear side 
conditions can be defined. 
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5 Generalized Stackel form for the Helmholtz equation 

Now we perform an analogous construction of eigenfunctions for a Helmholtz-like operator, 
using the same generalized Stackel matrix S. We want eigenfunctions of the separated form 

N 

* = II ^^\xi). We take separation equations in the form 

+ ftffidt* + \ vt(x l ) - ^s to (a; £ )A a ^ = 0, £ = l,...,N. 
Then we have the eigenvalue equations 

A' 



LP* = Y, T ^( d t + hd t + v e )V = XpV, = 1,...,N-1, 



and the side condition 

N 



L N V = J2 Tm ( d " + f& + = °- 



=i 



We take L 1 = H, X\ = E, so —\H is the standard Hamiltonian operator. 
Let 

X e = dj + f e d e + v t , Y t = 8j + f e 8 e . 
We need to compute the commutator [L a , U] for a = 1, . . . , N — 1, j = 1, . . . , JV. Now 

V i J \ k J i,k 

i,k i,k 

VL"=(YPk j) H k 2 Xk) (Ep?^- 2 ^) =Ep^P^ H ^X kXi 

V k J V i J i,k 

+ Epk )H k 2 MP < r ) Hr 2 )Xi + 2Yp k 3) Hfd k { P i t ) H^) dkXl , 



so 



[L a ,L j ] = £ (prWV) - pfy^prHf)) *fx k 

i,k 

i,k 



+ 2 Y(pi di (Pk H k 2 ) ~ Pi d i{pk H k 2 )) H^diXk- 



Using (4.7) we can establish the identities 

N 

_ AN) IT- 2TT- 2 ( JO) (j) (a) TT-lirs n \f> 

- Pk H k H i l^yPh Pi ~ Pi Ph ) H h \Vi a h)Oi 
h=l 



10 



P. Broadbridge, CM. Chanu and W. Miller Jr. 



N 

for any function F, and 

JV 

_n n (N)jT- 2q rr-2 /„(«) Jj) n 0) -("h fj— 2n _ 
- 2 Pi ^fc P i ~ P h Pi ) H h d i a h 

h=l 

+ P^H^Hf (pf > - - P^pf)^^ . 



Thus, 



AT / / N 



i=l V V h=l / 

x (-i e - inn?** + + (e or - 

x (E - Pfpf )^ 2 a t a h X\ L N = F ai L», 

where F a j is a first-order partial differential operator. 

We see that there is no obstruction to lifting our classical nonregular separation to the 
operator case. A difficulty occurs, however, when we try to write the pure operator part of H 
as a Laplace-Beltrami operator on a Riemannian manifold. Then there is an obstruction, a gene- 
ralized Robertson condition, to be worked out. Also we need to examine the effect of permitting 
i?-separation. 



6 Maximal nonregular separation as regular separation 
with a side condition 

Another way to approach the classical Hamilton-Jacobi problem is to use the Kalnins-Miller 
method for variable separation [17, 25] and consider maximal nonregular separation as regular 
separation with a side condition. (Here the nonregular separation is maximal in the sense that 
with a single side condition the number of separation constants is the maximum possible for 
nonregular separation, i.e., just 1 less than that for regular separation.) We look for additively 
separable solutions of the equation 

N 

Y J H~ 2 u 2 l +V = E (6.1) 

i=l 

with the side condition 



N 

Y,L- 2 u 2 + W = 0, 

i=i 



(6.2) 
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i.e., solutions u such that Uij = for i ^ j. From (6.1) we find 



■it 



'33 o Tj-2 ' 

2H j Uj 
and from (6.2) 

2LJ\ ■ 

These expressions must be equal modulo the side condition (6.2), so 

v, ± E, asv _ gj + E, ^ r y / E v \ 

for some functions z/j. Similarly, equations derived from Uf/fc = for j ^ k must hold modulo 
the side condition. Requiring that all of the above equations hold identically, i.e., requiring that 
we have regular separation modulo the side condition we eventually obtain the conditions that 

1. There are functions uiji = uiij for all j ^ i such that 

-2 



-djlLf + & ^ r d £ Lf + &e yi r d j Lf + u j£ Lf = (6.4) 

and 



J 2 1 J ' 



8.LJ 2 d e Lf 

-W jt + Wi- 1 ^- + Wj— + = (6.5) 
L e L- 

for all i = 1, . . . , N. (If all the Li are nonzero, this means that the L~ 2 are in confor- 
mal Stackel form, i.e., an arbitrary function times a Stackel form matrix, and that W is 
a conformal Stackel form potential [23, 24].) 

2. There are functions Tj, j = 1, . . . , N such that 

S ! = ^ ! +rA- 2 , (6-6) 
for all i = 1, . . . , N, and 

^ = S + ^- (6-7) 

3 3 

3. Let Cij be the second-order differential operator acting on functions / by 

d,H~ 2 di\nHT 2 
CaU) = dijf ~ -i^dj - —jpf-djf. (6.8) 

There are functions = fi£j for all j 7^ £ such that 

C je {Hr 2 )=» jt Lr 2 , (6.9) 

and 

C jt (V) = HjiW (6.10) 
for all j,£ = l,...,N, with j / 1. 
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7 Maximal nonregular separation =>- generalized Stackel form 

We have shown that the first row of the inverse of a generalized Stackel matrix with arbitrary 
iV-th column is an orthogonal metric H^ 2 whose associated Hamilton-Jacobi equation admits 
nonregular separation on the hypersurfaces given by the level set C N = of a function C N , 
quadratic in momentum variables, which is a constrained first integral on the same level set and 
whose components are the iV-th row of the inverse of the generalized Stackel matrix. Now we 
prove the converse, i.e., that if an orthogonal geodesic Hamiltonian is separable in orthogonal 
coordinates on the level set of a quadratic first integral, then it is a row of the inverse of a gene- 
ralized metric and the quadratic coefficients of C N are the iV-th row of the generalized Stackel 
matrix. Our starting point here is the geometrical framework of regular separation of variables: 
a complete separated solution of the Hamilton-Jacobi equation is a foliation parametrized by iV 
parameters for the integral manifold of the distribution generated by the N vector fields D{ = 
d x i + Rid Vi where the Ri are determined by the condition that the D{ are tangent to H = const. 
The classical Levi-Civita conditions [22] represent the integrability conditions of the distribution. 

In our case we need this distribution to be integrable only on the submanifold S defined by 
C N = 0. We also need the vector fields to be tangent to the submanifold S (closely related to 
the compatibility of the side condition). Our first step will be to write the differential conditions 
that in this case play the role of the Levi-Civita condition for regular separation. In this case 
they mix H and CJ^ . Moreover, they include also the condition that D{ are tangent to C N = 0. 
We will show that these equations are exactly the equations for nonregular separation derived in 
Section 6. Then we will relate these conditions to the existence of some more intrinsic geometrical 
object for the eigenvalues of Killing tensors (associated with quadratic in the momenta constants 
of motion), whose integrability conditions are equivalent to the Levi-Civita conditions. 

N 

Finally, we will construct the family of quadratic "first integrals" C h = ^2 T h ^u 2 commuting 

3=1 

and constant for the motion on C N = diagonalized in the coordinates we are considering, and 
show that the inverse of the matrix of the components T h i is a generalized Stackel matrix. Then 
we will extend our analysis to prove the corresponding results for multiplicative separation or 
i?-separation of the Helmholtz or time independent Schrodinger equations. 



7.1 Differential conditions for nonregular separation 
on a quadratic first integral leaf 

As in Section 6 we consider a natural Hamiltonian in orthogonal coordinates x = (x l ) on the 
cotangent bundle of a iV-dimensional Riemannian manifold Q with Hamiltonian 

N 

n = £ 1 = Y,Hr 2 p 2 + V(x), (7.1) 

and a function C N which is also quadratic in the momenta (pi) and diagonalized in the same 
coordinates 

N 

8=1 

where pf 1 are the eigenvalues with respect to the metric H~ 2 , i.e., = H~ 2 . We want to 
study the existence of separated solutions u of the Hamilton-Jacobi equation 

N 

H~ 2 u 2 + V(x) = E, Ui = d iU , (7.2) 

i=i 
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with the side condition, or constraint, 

N 

r' X! > 

=1 



2>fVii? + W(x) = 0. 



Recall that the conditions for nonregular separation are (6.4), (6.5), (6.6), (6.7), (6.9), (6.10) 
for Li = pi H^~ 2 and some functions u>£j, Tj, Pji- However, in this case, we can avoid the 
introduction of additional unknown functions, since we can solve the equation C N = with 
respect to a momentum variable Ui (as for instance u\ = u\{x^ ,u a ) with a = 2,...,N) so 
an expression vanishes on C N = if and only if it vanishes for all u a for a = 2, . . . ,N after 
the substitution of u% by u\{x 3 ,u a ). This simplifies very much the task of finding equivalent 
conditions such as the link with generalized Stackel matrices and it is possible only because we 
are assuming orthogonal coordinates. 

The conditions (6.4), (6.5), (6.6), (6.7), (6.9), (6.10) are equivalent (supposing without loss 
of generality that p^' 7^ 0) to imposing that for 

w N J^H-z 

v 2 - - V Pa 01 v 2 (7 -1) 

Pi H \ a=2 Pi H l 

the expressions (6.4)-(6.10) vanish for all values of u 2 a . Even easier, inserting (7.3) in (6.3) we 
find that the coefficient of v 3 vanishes and equating coefficients of u 2 a we get 



WlT-2 fi.TT-2 fl.„Wff-2 » . TT~2 fl.U/ 

OjV. (7.4) 



d jf g>H? P r j d jH - 2 pT'd^ 2 _d 3 w p 3 



n^U- 2 „M m 2 n^H- 2 H7 2 W W 

p a H a p a n a p x H 1 Pi n \ 

Similarly we find 

q,(H- 2 ) _ C tJ {H^ 2 ) _ C l3 (V) 
p^Ha 2 p[ N) H^ 2 " W ■ 



(7.5) 



In all these equations we follow the convention that the vanishing of a factor p a N ^ or W in 
a denominator in one of expressions (7.4), (7.5), implies that the numerator vanishes. 

Note that (7.5) and (7.4) give conditions for separation of the geodesic Hamiltonian with 
V = W = 0, with additional conditions that V and W must satisfy. Supposing for simplicity 
for the moment V = W = 0, we can rewrite (7.5) and (7.4) as 

(7.6) 



C ij { H a 2 ) 


AN) 
Pa 


(H 1 


2 ) 


H a 2 


(N) 
Pi 


Hi 2 




AN) 

d pa - 


( (AO 
( Pi 


pi N) \ 


diH~ 2 


1 (AO ~ 

Pi 


I (A r ) 
\Pl 


p{ N) ) 


H a 2 



AN) ( (N) \ -j rr-2 

+ E^- [ 1 - Pi— ) 9iH i (7 7) 



which are the necessary and sufficient conditions. Equations (7.7) can be interpreted as a 

(JV) 

first-order system in the N — 1 unknowns B a = ^fm- 

Pi 

Proposition 1. A geodesic Hamiltonian % admits nonregular separation on the submanifold 
C N = in a giv< 
1, . . . , N) satisfy 



(JV) 

C N = in a given orthogonal coordinate system if and only if the functions B^ = %vt (^ 

Px 



3^ = (Bi - Bu) 8 ^- + B k (l - Bi) 9 ^. (7.8) 
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Proof. By rewriting the necessary and sufficient conditions (7.6) and (7.7) in terms of the N 
functions B k = B\, B a with B\ = 1 we get 

°ij ( H k 2 ) R °ij i H l 2 ) f>-7 q\ 

— 2 — =H k 2 (7.9) 

and (7.8). However, it is a straightforward calculation that (7.9) is a differential consequence 
of (7.8), indeed we have 

didyB* - djdiBk = (Bi - Bj) \ Clj{ L_ k 2 ^ - B k Cij< f_l A . M (7.10) 

Note that the system of equations (7.8) coincides with (4.9), verifying again that separation 
with a generalized Stackel matrix is nonregular. 

Proposition 2. A natural Hamiltonian % admits nonregular separation on the submanifold 
C N = in a given orthogonal coordinate system only if the ratios 

Cij(H~ 2 )H^ 2 
Cij(H 1 2 )H a 2 

are independent of i and j and the eigenvalues of the quadratic function C are proportional to 
them 

_pi N) _ C %J (H- 2 )H^ 

p[ N) " C tJ {H^)H^ (7 - iij 

Remark 3. The function C N is naturally defined up to a multiplicative factor / on Q: if an 
expression is zero on C N = then it is also zero on fjC = for all functions on T*Q, but, since 
we are interested on quadratic in the momenta functions, we can normalize / on Q. (From the 
Stackel matrix point of view this corresponds to the multiplication of the iV-th column by /. 
Hence, equations (7.9) determine the unique (up to a factor) quadratic hypersurface where 
separation could occur. Indeed, equations (7.10) are the complete integrability conditions for 
the first-order PDE system (7.8). These conditions are identically satisfied for all {Bk,x l ) on 
an open subset of C 2n only if dj(H^ ) = 0, that is only if regular separation occurs. However, 
a single solution B k = Bk{x % ) could exist, provided it satisfies the original conditions (7.8), 
that is if it takes the form (7.11), or with all B^ equal, that is for C N = 1-1. In this case we 
get separation of the null equation T~L = (in which case a column of the Stackel matrix was 
arbitrary, so the result is consistent with the known one). 

Let us examine the relation of this kind of separation with the generalized Stackel matrix. 

Theorem 1. Suppose the natural Hamiltonian (7.1) admits N — 2 other functions C 2 , . . . , C 1 ^' 1 , 
quadratic in the momenta such that 

1) C 1 = H, C 2 , . . . , C 1 ^^ 1 , C are pointwise independent, 

2) C 1 = H, C 2 , . . . , £ N ~ 1 ,£ N are constants of the motion, modulo C N , that is 

{H,C k }\ cN=0 = 0, 

3) the quadratic terms of C 1 = 7~L, £ 2 , . . . ,£ Ar_1 , C N are diagonal in the coordinates (x % ), that 



■is 

N 



^ = £/f^ 2 K 2 + ^(x)), P v = i v, 
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Then, 



1) djVi = for j 7^ i and the eigenvalues pf of satisfy the following generalization of the 
Eisenhart conditions 



2) i/ie inverse matrix ofT h = (p^pH i 2 ) is a generalized Stackel matrix with the last column 
made of arbitrary functions of N variables ai(x), . . . , a n (x); 

3) the Hamilton Jacobi equation (7.2) admits a maximal nonregular separated solution on 
C N = depending on N — 1 parameters; 

4) the N functions C 1 = %, C 2 , . . . , C N ~ l , C N are in involution on S, that is 

{C l ,C k }\ C N =Q . 

Proof. 1) By inserting (7.3) in the Poisson brackets 



_ / N 

i \k=l 

N \ 

k=l I 



we get that diV k = for i ^ k and that p^p satisfies the generalized Eisenhart conditions (7.12). 

2) Since the quadratic diagonal functions are pointwise independent, the matrix T" has 
a nonzero determinant and therefore it admits an inverse s k h such that 

T ik s kh = 6{ (7.13) 

(we partially follow the proof of equation (4.7) and a proof of the Eisenhart theorem given 
in [11]). We want to show that <9;s^ = for all i ^ h and £ / N if and only if pf ] = T £j /T lj 
satisfy the generalized Eisenhart conditions (7.12). Differentiating (7.13) with respect to x l we 
get 

(^T ik s kh + T ek d lSkh ) = 0, 

k 

that is 

{{Hfdipf + p k e) diH- 2 )s kh + p k e) H- 2 diS kh ) =0, £^1, (7.14) 

k 

{^H k 2 s kh + H k 2 d lSkh ) =0, l = \. (7.15) 

k 

Adding (7.15) multiplied by —pp to equation (7.14) we get 

E m 2 *P? + {Pf - Pf)d^ 2 )s kh + {pf - pf)H~ 2 d iSkh ) =0, t + 1. (7.16) 
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Let us suppose that diSkh = for all i ^ k and h ^ N and diS^N = 9i0<hi by multiplying (7.16) 
by T h \ we get 



Hj 2 d iP f + { P f - P^mf + E ((^ " P?) H k 2 dia k ) = 0. (7.17) 



k h 

that is 

A' 



Writing (7.17) for j = 1 and j ' = a and by combining them in order to eliminate the common 
factor Ylk {{pf ~ Pi jH^didk) , we get (7.12). Conversely, let us suppose that (7.12) hold; this 
means that the quantities 

H k d *Pk + [Pk ~ Pk ) d i H k 



J^Nk 

are independent of k, that is H k 2 dip k ^ + (p k ~ pf)®iH k 2 = T Nk Q[. Therefore (7.16) becomes 

T Nk QUk h + E {{Pf ~ P?)H k 2 d lSkh ) = 5M + E ((P? " P?)^ 2 ^) = 0. 
fe it 

If we fix the values of h and i, then the above equations are a linear system of N — 1 equations 
(£ = 2, . . . ,N) in iV — 1 unknown diS^h k = 1,...,N, k ^ h). For h ^ N the system is 
homogeneous and the matrix of coefficient M = yp k — p\ )H k has nonzero determinant 
(see [11]). Hence, the only solution is diSkh = for h ^ N. 

3) If equations (7.12) hold, they hold in particular for £ = N and we get (7.7) which are 
proven to be equivalent to conditions (7.8) that are necessary and sufficient for the nonregular 
separation with the side condition C N = 0. Further, the conditions (7.4) on V, W are exactly 
the consistency conditions that must be satisfied by V, W in order that {T-L,£, N } = hold on 
C N = 0, and conditions (7.5) are exactly the integrability conditions for V, W that must be 
satisfied if {H, C N } = is to hold. 

4) This is essentially the computation in Section 4. ■ 

Theorem 2. // the N pointwise independent functions 

£} = U,C 2 ,...,C N ~\C N 

are in involution on C N = and admit common eigenvectors, i.e., if the quadratic forms in 
the momenta associated with — p^'T-L have N common eigenvectors, then the eigenvectors are 
normal: there exist orthogonal coordinates (x l ) such that the N functions are simultaneously 
diagonalized. 

Remark 4. This is essentially the converse of the results of Section 4. It shows that separation 
with a side condition implies the existence of a generalized Stackel matrix that defines the 
separation. 

Proof. Apply Theorem 7.5 of [2] with M a j, = K , where K N is the symmetric 2-tensor asso- 
ciated with C N . ■ 

Remark 5. The functions are defined up to multiples of C N and the conformal Killing 
tensors are defined up to multiple of the metric tensor. Indeed, C e + f e C N have the same 
eigenvalues of C e (provided that the £ e have common eigenvectors) and they are in involution 
when restricted to C N = 0: 

{C" + fC N ,C k + f k C N } = QC N . 
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Now we close the logical loop. Suppose the conditions (7.8)-(7.10) of maximal nonregular 
separation are satisfied and consider the linear system of equations 

diPj + {pj - pi) % H ^ 2 = Bj i^dipx + {pi - pi) % ^} 2 j , i, j = 1, . . . ,N, (7.18) 

for the unknowns (pi, . . . , pjv)- Remark 5 suggests that this system has solutions for any choice 
of p\. Therefore, we choose the function p\ arbitrarily and consider (7.18) as a system of 
N(N — 1) independent differential equations for the unknowns (p2, ■ ■ • , Pn)- It is then a straight- 
forward exercise to verify that the integrability conditions for (7.18) are satisfied identically, due 
to the nonregular separation conditions. We already know the two solutions (1,1,..., 1) and 
(i?2, • • • , Bjst). Indeed about any regular point xo we can find a unique solution such that 
Pj(xo) = f3j, 2 < j < N for any choice of constants f3j. For example we could take p± = 1 
and define N — 1 linearly independent solutions by choosing the f3j accordingly. Then we could 
take p\ = and find another solution such that the full N solutions form a basis. This basis 
determines a generalized Stackel matrix. 

It remains to consider the potential terms, the given functions V and W. Using the func- 
tions pi of the generalized Stackel matrix computed above we define functions 

N 

£ £ = E^rV + ^ w , £=i,...,n-i, 

i=l 

where V = but the remaining will be determined by imposing the requirement that 

there exist functions Fj linear in the momenta and such that 

{%, C N ] = F N £ N , {H,C e } = F e C N . 

By working out the Poisson brackets, it is straightforward to verify that equations (7.4) are 
exactly the consistency conditions for the existence of Fn and (7.5) are exactly the consistency 
conditions for the Fg and the integrability conditions for the existence of the potentials 

Now note that since the generalized Stackel matrix is invertible we can uniquely determine 
functions Wj(x) such that 

N 

W W=J2viP?Hr 2 , j = l,-..,N. 

8=1 

With this observation we have verified all of the assumptions of Theorem 1. This proves the 
following result. 

Theorem 3. // a natural Hamiltonian % admits maximal nonregular separation on the sub- 
manifold C N = in a given orthogonal coordinate system, then the system is separable with 
a side condition and there exists a generalized Stackel matrix for the separation. 

All of these results for additive nonregular separation of Hamilton-Jacobi equations extend to 
multiplicative nonregular separation for the Helmholtz (or Schrodinger) equations on the same 
manifold, but with an obstruction. If a system is nonregular i?-separable for the Helmholtz 
equation then it is nonregular separable for the Hamilton-Jacobi equation. However, if the 
Hamilton-Jacobi equation admits nonregular separation in some coordinate system then there is 
a "generalized Robertson condition" to be solved to determine for which potentials the Helmholtz 
equation admits i?-separation in these coordinates. On the other hand, we can always find some 
family of vector potentials for which the Helmholtz equation does admit nonregular i?-separation 
in the coordinates. That is, it is easy to show that given a nonregular separable system for 
the Hamilton-Jacobi equation, we can construct a family of vector potentials for which the 
corresponding Helmholtz equation is nonregular i?-separable. We give examples in Section 8. 
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8 Examples of nonregular separability 

for Hamilton— J acobi and Helmholtz equations 

8.1 Examples of restricted regular separation 

The simplest examples of Theorems 2 and 3, and their extensions to the Helmholtz equation, are 
those for which regular orthogonal separation already occurs. Thus the "generalized" Stackel 
matrix (4.1) is a true Stackel matrix. Consider first the separable Hamilton-J acobi equation in 
the form 

N 

1=1 

Here, the metric in the orthogonal separable coordinates x 1 is 

N 



ds 2 = ^2H?{d 



i=l 



There is an associated Stackel matrix and constants of the motion C^, 2 < t < N, where % = C 1 . 
We apply the side condition C = C N — Xn = where A at is a constant scalar. Note that C 
is also a constant of the motion. Now the separation equations become 

N-l 

u 2 i + (v i (x i )+X N s iN (x i )) + Y J s ij {x i )X j = 0, i = l,...,N, X x = -E. 

3=1 

Our "restricted constants of the motion" become 

N 

Ct'^J^^itf+Vj + XNSjN), 
3=1 

so we have modified the potential. The side conditions are satisfied automatically. Effectively, 
we have restricted our Hamiltonian system to the N — 1 dimensional hypersurface C N = Aat and 
seen that the result is again a Hamiltonian system. Under some circumstances this has a simple 
interpretation as a separable system on a manifold of one less dimension. 

Example 1. Consider the Kepler Hamiltonian 

U = C 1 =p 2 x +p 2 y +p 2 z + r = yV 2 + y 2 + z 2 

in real Euclidean space. The system is separable in spherical coordinates r, 9, <j> where 

x = r sin 9 cos <j>, y = r sin 9 sin </>, z = r cos 9, 
r > 0, < 9 < 7T , < < 2tt. 

Here 

r 2 rv r 3 

r z r sm z 9 v 

We choose the side condition £ 3 = As- Then the reduced Hamiltonian becomes 



r z r z sin 9 r 
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In terms of the new variables X = r cos 9, Y = r sin 9 this becomes 

nil 2,2 a , ^3 

n =Px+Py + VWT^ + y~^ 

a regular separable system in the first quadrant of the Euclidean plane. In the special case A3 = 
the singularity on the X-axis disappears and we can extend the system to the full punctured 
plane. Note that this construction is conceptually distinct from simply restricting a trajectory 
confined to a plane. 

A similar construction works for the Helmholtz equation. We want eigenfunctions of the 

N 

separated form \I> = Y\ ^^(x^). We take separation equations in the form 
i=i 

dfo + fe{x e )d^ + (^v e (x e ) - s m {x e )X N - J2 ^{x^X^j * = 0, £ = 1,...,N. 
Then we have the restricted eigenvalue equations 

N 

= T ^i d t + Wt + v e~ sw\n)* = A^*, p = l,...,N-l, 

e=i 

and the side condition 

N 



We take L 1 ' = H' . Effectively, we have restricted our separable quantum system to an eigenspace 
j* : L N ^> = An*} of a symmetry operator. The restriction is again a Hamiltonian system. 
Again, under some circumstances this has a simple interpretation as a separable quantum system 
on a manifold of one less dimension. 

Example 2. Consider the hydrogen atom Hamiltonian, a constant multiple of 

H = L 1 = d 2 x + d 2 + d 2 + r = vV + y 2 + z 2 

in real Euclidean space. The system is separable in spherical coordinates r, 0, (p: 

H = d 2 r + -d r + \l 2 + -, L 2 = dl + cot 9d e + -r^L 3 , L 3 = dl. 

r r z r sin v 

We choose the side condition L 3 * = A3*. Then the reduced Hamiltonian becomes 

H =d r + -d r + -^d 9 + -^-d e + . 

r r z r z sin # r 

Now we make an ^-transformation * = R@ where R = — r=. Then the eigenvalue equation 

rv sin 

H'ty = £* becomes 

(d 2 + \d 2 + ^i- + i + ^]@ = Ee, d 2 e = x 3 @. 
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Now write = <&(</>) E(r, 6) . In terms of the new variables X = r cos 8,Y = r sin 8 the energy 
equation is 

a regular separable system in the upper half Euclidean plane. If A3 = —1/4 this equation can 
be extended to the punctured plane and regarded as a modification of the hydrogen atom in 
the plane. Note, however, that due to the R factor, the usual Li normalization for eigenstates 
in Euclidean space doesn't restrict to the usual L2 normalization for eigenstates in the plane. 
Thus, in analogy with the virial theorem for the mapping between the Coulomb and pseudo- 
Coulomb problems [36], one needs to check that the bound state spectra are preserved under 
the restriction. 



8.2 Nonregular separation in 2D and a "no go" theorem 

Now we take up the issue of true nonregular separation in a coordinate system in which regular 
separation is impossible. Let us first look at 2D examples, a very special case. By using the 
facts that we can always replace a separable coordinate by an invertible function of itself, and 
we can perform linear transformations on the separation constants without changing the system, 
we can always assume that the generalized Stackel matrix looks like 

q _ (1 A(u,v) 

We must require that A — B 7^ so that the matrix is invertible and that AB 7^ so that the 
metric is nondegenerate. Further, by making use of the equivalence relation (4.10) we can put 
the matrix in canonical form 

< i 1 1 \ -ti, • r W / -1 

d = I, , , with inverse 1 — ' 



1 f(u,v))> /-iV" 1 1 

where / = B/A. This form will lead to true nonregular separation unless / can be factored 
as f(u,v) = U(u)V(v), because in that case we can take B = V(v), A = 1/U(u) and find an 
equivalent true Stackel matrix. It is easy to see that this form leads to true nonregular separation 
for the Hamilton-Jacobi equation. However, in the 2D case the construction fails for Helmholtz 
equations. 

Theorem 4. For a 2D manifold the Helmholtz equation A^/ = A^ never admits true nonregular 
separation. 

Proof. In the 2D case, the Laplacian takes the form 

/ (q + IU 8 \_J_( 8 _l_f ld 

u uu \ n j, U U I j> -1 I '-'vv n j, u v 



/ — 1 y 2 / y /-iv 2 / 

For nonregular separation it is necessary that f u / f is a function of u alone and f v / f is a function 
of v alone. Thus / = U(u)V(v) and the system admits regular separation. Similarly if \P = e R ® 
for some fixed R, then i?-separation implies that there exist functions U(u), V(v) such that 

R u = -\ f f + U(u), R v = ^ + V(v). 

Since d v R u = d u R v it follows easily that d uv In / = 0. Thus / = JJ(u)V(v) and the system must 
admit regular i?-separation. ■ 
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However, true nonregular separation may occur when vector potentials are included, even 
in the 2D case. The self-adjoint Schrodinger equation for a non-relativistic charged particle 
influenced by a classical electromagnetic 4-potential, still does not admit true nonregular sepa- 
ration. Nonetheless true nonregular separation does occur when the vector potential gives rise 
to skew-adjoint first-order terms, as in a diffusion-convection equation for transport of solute. 

8.3 Quantum particle in IR 2 with a magnetic field 

In the Schrodinger picture with Cartesian coordinates, the eigenvalue equation for the quantized 
Hamiltonian of a non-relativistic charged particle interacting with a classical electromagnetic 
field, is (e.g. [4]) 



i_L V 2 ^ - y [A.V + V.A(x)]* + *(x)tt = A*. 



(Ai,A 2 ,0) is the vector potential for the magnetic field and ^ is the scalar potential for the 
electric field. We consider the equivalent eigenvalue problem 



= + <9f* + - [a l {yi)d 1 + a 2 (x)<9 2 + V.a(x)] ^ + $(x)^ = EV. ( 
We extend this equation to a general curvilinear orthogonal coordinate system. As before, 



U) 



H^ 2 dl + H 2 2 d 2 2 + #r 2 di log - Hfd 2 log ^d 2 



In a general orthogonal coordinate system, the scalar divergence, involving Christoffel symbol T, 
of vector a, is 



Ba m • 1 8 

diva=^- + a m n rr " 



dx m im y /gdx i 

did* + a 1 [dx log Hi + d x log H 2 ] + a 2 [d 2 log Hi + 8 2 log H 2 ] . 



The eigenvalue equation for the Hamiltonian is 



= tff 2 <9 2 * + H^d'fo + 



2q2, 



#i dilog 



El 



Ho 2 d 2 \og 



El 



d 2 V 



+ - [^(x)^ + flia^x) + a 2 (x)<9 2 + <9 2 a 2 (x)]l' 



+ - [a}d! log(ifi) + a 1 8 1 log(H 2 ) + a 2 8 2 log(#i) + a 2 8 2 log(tf 2 )]* + $(x)* = EV. 
Consider i?-separation 

e = *e- R W = Q^{x r )Q^{x 2 ). 

Then 



H^df + H^d 2 + H^ { di log ~ + 2d x R \ d x + H^ 1 { -8 2 log ~ + 2d 2 R j> d- 



H 2 



d 2 R + (dxR) 2 + {d 1 R)d 1 log 

3\R + (d 2 R) 2 - {d 2 R)d x log 



H 2 



Hi 



e 



El 

Hi 

El 

Hi 



e 

e + $(x)e 
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+ l - [a l dxQ + a 2 d 2 Q] + % - [c^dxR + a 2 d 2 R + fra 1 + d 2 a 2 ] 9 
+ V[di logfli + dx \ogH 2 ]Q + ^a 2 [a 2 log^i + <9 2 log# 2 ]9 = EQ. 
Assume the separation equations 

d 2 Q + [hix 1 ) + igiix^dxQ + [Mx 1 ) + ^(x 1 )] - Suix^E] 9 = 0, 
d 2 2 + [f 2 (x 2 ) + ig 2 (x 2 )]d 2 Q + \[v 2 {x 2 ) + iw 2 (x 2 )} - S 21 (x 2 )E] 9 = 0, => 
T u dfe + T 12 c> 2 2 9 + T n [h + %]ai9 + T 12 [/ 2 + ig 2 )d 2 Q 

+ [Th[di + iwi] + T\ 2 \v 2 + w 2 ]] 9 = EQ, 
T 21 d\Q + T 22 a 2 2 9 + T 2l [h + %]9i9 + T 22 [/ 2 + i 52 ]d 2 9 

+ p2i[t>i + iwi] + T 22 [u 2 + iw 2 ]] 9 = 0, 

where matrix T is the inverse of generalized Stackel matrix S. In order for the above eigenvalue 
equation to be identified with the Schrodinger equation, we again require T±\ = H^ 2 and 



T\ 2 = H 2 2 , and in addition, 

j^x 1 ) = di log ^ + 2diR, f 2 {x 2 ) = -8 2 log ~ + 2d 2 «, 
Tn«i + T 12 t; 2 = $ + H{ 2 8fR + fl^Sgfl + H^idxR) 2 + H 2 2 {3 2 R) 

+ H^ 2 {d 1 R)d l log ^ - i/ 2 - 2 (a 2J R)a 2 i og ^, 



gi[x 



-ti 1 a , 



92 (x 2 ) 



1 



:-W 2 a > 



(8.2) 



Tnwi + Ti 2 ?i; 2 = - [a l d x R + a 2 d 2 R + ^a 1 + d 2 a 2 ] 

+ V [9! log ff! + «9i log ff 2 ] + X -a 2 [8 2 log ffi + 8 2 log fr 2 ] . (8.3) 

From (8.2) we can easily deduce that d\d 2 R = and d\d 2 \og(H 2 / H\) = 0. Hence 

H 2 /Hi =ni(x 1 )n 2 (x 2 ) and R = R 1 (x 1 ) + R 2 (x 2 ) (8.4) 
for some functions III, IT 2 , R\, R 2 . Now assuming the canonical form with H\ = 1, we have 



a 1 =2g 1 (x 1 ) and a 2 = 2g 2 (x 2 )U 2 (x 2 ) II^x 1 ) , 

whence it follows by separation of variables in (8.3), that 

Uj[wi - gidiRi - digi - gidi loglli] = 74, 
II2 2 [w 2 - g 2 d 2 R 2 - d 2 g 2 + g 2 d 2 log n 2 ] = -74, 

where 74 is constant. By integration, 



i l = 2nrV Kl 



75 + j e Rl [wtUx -^U^]dx 

76 + J e R2 [w^ 1 + 74 n 2 ] dx 2 



1 2 = 2n 1 - 2 n 2 - 1 e -- R2 



with 74, 75 and 76 constant. As a standard test, we first recover the known possibility of 
separation in polar coordinates in E 2 , with 



x 2 = e, 



n 1 = r, 



n, 



1. 
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In (8.2) we may choose f\ = f 2 = without any loss of generality to the class of allowable scalar 
and vector potentials. Consequently, from (8.2), R\ = — \\ogr and R 2 = 0. This yields the 
potential fields that admit separation of the Schrodinger equation 

a 1 = ai(r) = 7 5 r~ 1/2 + 2 74 r _1 + r~ 1/2 J r l/2 wi(r)dr, 

a 2 = a 2 (6)r- 2 = 2 76 r~ 2 + 2r~ 2 J 74 + w 2 (6)d0, $ = vi(r) + r~ 2 v 2 {6) - ^r -2 . 

The above form for the scalar potential agrees with that known to allow separation of polar 
coordinates in the Helmholtz operator [19]. When the vector potential a in (8.1) takes values 
in M 2 , the separation condition (8.4) again leads to the conclusion that in this case, generalized 
Stackel matrices do not lead to new separable systems. 



8.4 Solute transport 



The possibilities of new separable equations are broadened if we replace the pure imaginary 
coefficients of first-order terms in (8.1) by first-order terms with real coefficients. Consider 
a solute diffusion-convection equation in Cartesian coordinates, 

3 t y = = df^ + di^ - [V(x)di + q 2 (x)d 2 }^ - /x(x)«. 

In this application, ^ represents the solute concentration, q(x) represents a steady velocity 
field of the solvent, and /u(x) is an adsorption coefficient for removal of solute by the solid 
substrate of a porous medium or by another (usually solid) component of the mixture. If ^(x) 
satisfies H^> = E^>, then e £ ''^ / (x) is a solution of the time dependent solute equation. In 
standard applications, E < 0, unless the solid component of the mixture is releasing solute 
(/i(x) < 0) rather than adsorbing solute. After making the replacements q = (— i/2)a and 
//(x) = — <5(x) — (i/2)V.a, the solute equation is directly analogous to the Schrodinger equation 
with magnetic field, except that now q 3 must be real, corresponding to a 3 being pure imaginary, 



-iA>, A 3 



-2q 3 , 



d t i& = = d 2 ^ + flftf + - [A^^di + A 2 (x)<9 2 + V.A(x)] $ + $(x)* 
In a general orthogonal coordinate system, 

= H{ 2 d 2 ^ + H^ 2 d 2 ^ + 



tfr^iogf 2 

Hi 



H 2 2 d 2 lo g H2 



d 2 V 



+ Y - [A 1 ^ + ftA^x) + A 2 (x)d 2 + d 2 A 2 (x)] * + 1 [A 1 d 1 log(fli) 
+ A 1 di log(H 2 ) + A 2 d 2 log(fli) + A 2 d 2 log(H 2 )] * + $(x)* = EV, 

consider i?-separation 

G = Ve- R M = ^ 1 \x 1 )¥ 2 \x 2 ). 

Then 



H^ 2 d\ + H^ z d' 2 + 



2q2 



5, log ^ + 2d x R ) <h + II 



-d 2 log^ + 2d 2 R\ih 



e 



5 1 2 J R+(ai«) 2 + (9ii?)5i log 



f?2 



e 
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dlR + {d 2 Rf - (d 2 R)d l log 



H 2 



e + $(x)G 



+ X - [A x d x ® + A 2 d 2 Q] + X - [A l d\R + A 2 d 2 R + d x A l + d 2 A 2 ] Q 
+ l -A x [dx log H x + dx log F 2 ]9 + ^A 2 [d 2 log #1 + <9 2 log H 2 ]Q = EQ. 
Assume the separation equations 

dfe + f 1 (x 1 )d 1 e+ [v 1 (x 1 )-s 11 (x 1 )E]e = o, 

d 2 2 + f 2 (x 2 )d 2 e + [v 2 (x 2 ] - S 21 (x 2 )E] 9 = 0, =► 

Tndfe + r 12 a 2 2 e + T n hdi® + r 12 / 2 3 2 e + [t xiVi + t 12 v 2 \o = Ee, 
r 21 dfe + r 22 <9|e + TnhdxG + r 22 f 2 d 2 e + [T 2lVl + t 22 v 2 ]o = o, 

where matrix T is the inverse of generalized Stackel matrix S. In order for the above eigenvalue 
equation to be identified with the solute equation, we again require T\\ = H± , T\ 2 = H 2 . 
Since the first-order terms of the target solute equation no longer have imaginary coefficients, 
the identification of first-order terms now leads to more general possibilities 



H 2 q l = -h {x 1 ) + d x log + 2dtR 

til 



Hh 2 = ~h (x 2 ) - d 2 log ^ + 28 2 R. 

-Hi 



1.5) 



By differentiating throughout, this implies 

did 2 log § + 2d!d 2 R - d 2 (H 2 q l ) = 0, 
-Hi 

This system is equivalent to 

Ad^R - d 2 (#iV) - d^Hlq 2 ) = 0, 



Ho 



-d x d 2 log -f + 2d x d 2 R - di (H^ q 

-Hi 



2d 1 8 2 log I 2 - d 2 (Hfq 1 ) + d x (H 2 q 2 ) = 0. 



Similarly, identification of the zero-th order terms also leads to a more general condition 
T n vi + T 12 v 2 = -n + H^dfR + H 2 2 3 2 2 R + H^ 2 {d x R) 2 + H 2 2 {3 2 R) 2 



+ H^ 2 {d l R)d 1 log ^ - H 2 2 {d 2 R)8 2 log ^ - q l dxR - q 2 d 2 R. (8.6) 



^ - H 2 \d 2 R)d 2 lo g §- 
-Hi -Hi 

In principle, the compatible velocity field q is recovered from (8.5) after substituting these 
expressions in (8.6) which determines the function R. 

The condition H 2 jH\ = IIi(x 1 )II 2 (a; 2 ), which is the key condition for the possibility of 
replacement by a regular Stackel matrix construction, is no longer true in general but occurs 
only if 



d 2 (H!q 1 )-d 1 (Hiq 2 ) 



0. 



In a description with Cartesian coordinates for Euclidean space, this would be exactly the 
restriction that the solvent velocity field is irrotational; V x q = 0. Consider the canonical form 



S 



-i 



1 



/ -1 
-1 1 



/(x) - 1 
Then (8.6) may be written as, 

-fvi +v 2 = -fd\R + d\R + f(d x R) 2 - (d 2 Rf - fd 1 R8 1 log(-/) + 2fd 1 Rd 1 log(l - /) 
- hd x R + d 2 Rd 2 log(-/) - 2d 2 Rd 2 log(l - /) + f 2 d 2 R + (/ - 

In looking for examples of genuine nonregular separation, we must consider metrics for which 
H 2 jH\ does not separate in variables x 1 and x 2 . 
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Example 3. We consider an example in 2D Euclidean space. The metric is 

ds 2 = dx 2 + dy 2 = l^l^ldu 2 + (1 - f(u, v))dv 2 , 
f{u,v) 



f(u, v) = (u + v + ^(u + v) 2 - 4) . 



Here 

x = (u + v) cos(0 — it), y = (u + v) sin(c/> — it), 

u + v + A/(n + w) 2 - 4 



= - ^ u + « + y/(u + v) 2 -4) 



2 arctan 



and I it + f I > 2. The Helmholtz equation with vector and scalar potential takes the form 

A * " T-i G7 8 " " u(u) ) + r~i (tt 8 - - v{v) )) * - B *' 

The separation equations are 

(d uu + 17(«) - E)¥^ (u) = 0, (fiU + V(v) - E)*W („) = o, 

with * = *W(it)*( 2 )(i;). 

Now r 2 = x 2 + y 2 = (u + v) 2 , from which the polar coordinate may be taken to be r = u + v. 
Then 

x = r cos{4> — u), y = r sin(^ — u), 
implying that the polar angle coordinate is 
eft — u = 9 (mod 2ir). 

Also 

r + \Jr 2 — 4 \ 



- (r + \A" 2 - 4) - 2 



arctan 



Note that ±0 is an additive component of each of the variables u = <ft{r) — 6 and v = r — <ft{r) + 6. 
Therefore physically relevant solutions must be periodic with period 2tt in each of the variables u 
and v. Taking the simplest case w 1 (= U) = and / (= V) = 0, direct separation is possible 
with R = 0, and from (8.5) 



(ffi)V = ~(H 2 ) 2 q 2 = l/V(« + v) 2 -4 = [r 2 - 4] 



-0.5 



The velocity field is defined on the exterior to the circle of radius 2 centred at the origin. The 
squared magnitude of velocity is 

,n 2 „ 1 r„2\2 n „ _ tr2/l\2 . rr2/„2\2 1 



(^^.ej + (g 2 re 2 .e 2 = H 2 (q l ) + H 2 { q 2 )' 

After taking E = — ui 2 , separated solutions for solute concentration may be combined in a Fourier 
integral 

/•oo 

* = / A^e^^ 2 * cos(lj[u + 5(u)])cos(uj[v + e(oj)])du), 
Jo 

with amplitude function A(uj) and phase functions 5(u), s(u). 
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8.5 Nonregular separation in more than 2 dimensions 



In dimensions greater than two, both true nonregular separation and i?-separation occur even 
without an added vector potential. 

Example 4. This is a 3D Minkowski space example. The metric is 

Aw 2 



ds 2 = dt 2 - dx 2 -dy 2 = - 



(u 2 — v 2 ) 2 



(du 2 + dv 2 ) + dw 2 = H 2 du 2 + H 2 dv 2 + H^dw 2 . 



Here, 



t 



w 



2wuv 

9 9 ' 



and we require u, v real and w > 0, \u\ > \v\. 

The Helmholtz equation in 3D Minkowski space is 



A 3 9 = EG, 



where 



2\ 2 



A, 



(u 2 - v 2 ) 
Aw 2 



Aw 1 



8w 



We look for i?-separable solutions 

6 = e R V = e fl tfW(u)* (2) («)* (3) M> e R = ^. 
By direct calculation we can establish the operator identity 

e- R A 2 e R — E = R- 2 {d uu + A 2 ) + H~ 2 (d vv - A 2 ) + H~ 2 ( 8, 

Thus the generalized Stackel matrix and its inverse are 

(o 1 I {u 2 -v 2 ) 2 \ 
Aw 2 
{u 2 -v 2 f 



WW o 



E 



S = 



1 



Aw 2 

(u 2 -v 2 y 

Aw 2 



( {u 2 -v 2 f (u 2 -v 2 ) 2 \ 

Aw 2 Aw 2 
(u 2 - v 



J 



\ 



In terms of the ^ functions we have 



/ ( u 2_ v 2) 2 / , 

L 1 ^ = I - v . J (d uu + d vv ) +[d ww - - 



L ^ = - 



Aw 2 
(u 2 -v 2 f 



y (u 2 -v 2 y 

— -1 - — 

Aw 2 Aw 2 

1 1 0/ 



* = E^, 



Aw 2 

The separation equations are 



(dun + d m ) -d vv \y = -a 2 *, lH = (d uu + d vv ) * = o. 



(3 UU + A 2 )*« =0, - A 2 )* (2) = 0, f 0™ - ^ - #W 3 > = 0. 

In terms of the functions the final separated solution of A3© = EQ is 



w 
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Example 5. Consider the Helmholtz equation in 3D Euclidean space A30 = EQ. It is well 
known that this equation is regular separable in exactly eleven coordinate systems. Moreover 
true regular ^-separation does not occur for any constant curvature space [15]. This example, 
due originally to Sym [8, 30, 34, 35] and presented from our point of view, shows that nonregular 
i?-separation with two side conditions occurs in orthogonal coordinates distinct from the usual 
eleven. Consider Dupin-cyclidic coordinates u, v, w such that 

b 2 cos u cosh v + (c cosh v — a cos u) w b sin u(a cosh v — w) 



x = ; , y 



acoshv — ccosu ' a cosh v — c cos u 

bsinhv(w — ccosu) 
a cosh v — c cos u 

Here a, b, c are positive parameters, with b 2 = a 2 — c 2 , c < a, and u E [0, 2tt). For any real v we 
require c cos it < w < acoshv, so that x, y, z are real. Since 

ds 2 = dx 2 + dy 2 + dz 2 =( b(aC ° Sh " } ) 2 du 2 +( ^-"^ ) 2 dv 2 + dw 2 

\ a cosh v — c cos u J \a cosh v — c cos u J 

= H 2 du 2 + H\dv 2 + Hldw 2 , 
these coordinates are orthogonal. Thus 
A 3 = d xx + d yy + d zz 

= „ I „ \d u (H^H 2 H 3 d u ) + d v (HlH^H 3 d v ) + OUH^H^d^] . 

H\H 2 H 3 

We look for i?-separable solutions 

9 = e R V = e^^(u)^ 2 \v)^(w) t e R = 1 

y (w — c cos it) (a cosh v — w) 

By direct calculation we can establish the operator identity 

e- R A 3 e R — E = FT/ 2 (d uu + ^) + H 2~ 2 ~ E )- 

Thus the generalized Stackel matrix and its inverse can be taken as 

/0 1 
5= 1 

V ~ H l 2 ~ H 2 2 

In terms of the functions we have 

= (tff 2 (d uu + ^) + H 2 2 (dw ~ + d ^ ) * = m > 
L 2 * = (d uu + ^ * = 0, L 3 ^ = (d vv - * = 0. 
The separation equations are 

+ ^ = 0, f a w - ^ *( 2 ) = 0, (d ww - e)^ = 0, 

fW(u)$( 2 )(«)$( 3 )W 
— c cos it) (a cosh v — w) 

Here we have 

[L 2 ,L l ] = F 22 L 2 + F 23 L 3 , [L 3 , L 1 ] = F 32 L 2 + F 33 L 3 , [L 2 , L 3 ] = 0, 
for first-order differential operators Fij, so we have ^-separation with two side conditions. 







H 2 ~ 2 


H 









\ I 


1 
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Example 6. Again we consider the Helmholtz equation in 3D Euclidean space but now we 
choose coordinates u, v, w such that 



1 / 1 / N • 

x = (v + w) cos u, y = —j= (v + w) sin u, 
V 2 v 2 



z = —^(v — w). 
V2 



We take a scalar potential 
U{u) 



V(x,y,z) 



The metric is 



(v + w) 2 



+ V{v) + W(w). 



ds 2 = dx 2 + dy 2 + dz 2 = - [{v + w) 2 du 2 + dv 2 + dw 2 ] , 
2 



(f + w) 2 ' 



W !) — H W — l - 



Thus 



a 3 + f = a** + a ra + d zz + y 



i 



i?1 HoH; 



+ d w (H 1 H 2 H^ 1 d w )] + 



lJ32-n3 

I7(u) 



[d u (H^H 2 H 3 d u ) + 8 v (HlH^H 3 d v ) 
+ V(v) + W(w). 



(v + w) 2 

We look for i?-separable solutions 

G = e R y = e R yU(u)& 2 \v)^ 3 \w), e R = (2[x 2 + y 2 })^. 
By direct calculation we can establish the operator identity 



e~ R A 3 e R - E 



(v + w) 



:{d uu + U{u) + I) + 2 [ d vv + 



V(v) 



A 2 



W(w) E x 

+ ( <>.r,r I 2 2 • A 2 



Thus the generalized Stackel matrix and its inverse can be taken as 



^ 1 2- 

2 (v + w) 2 

-1 -2 



\0 1 

The separation equations are 



1 \ 



(2 2 2 \ 

(v + u>) 2 

0-1 -2 
\0 1 , 



(a uu + u(u) + = 0, (d w + ^ - a 2 ) ^ = 0, 



W(w) 



+ A 2 ] ^ (3) = 0, e = 



*(1)( U )*(2)( W )*(3) 



2 2 ' * ) yjv + w 

Example 7. As a generalization of the previous example, the Euclidean space metrics 

ds 2 = {Ui{u)w + U 2 {u)v + U 3 {u)) 2 du 2 + dv 2 + dw 2 

for arbitrary functions Uj(u) all lead to nonregular separation for Hamilton-Jacobi equations, 
and in the quantum case, to nonregular ^-separation for velocity dependent potentials. The 
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construction of separation equations and of the generalized Stackel matrix is standard. As an 
instance, if 

x = F\(u)v + F2{u)w , y = F${u)v + F±{u)w, z = F§(u)v + Fq(u)w, 

for 

7-1 2 7-1 ~~ C0S3 U 7^ 

F\ = sin u, F% = — — , F3 = sin u cos u, 



+ sin 2 n 



„ sin it (l + cos 2 u) „ ^ — sin 2 n 

F 4 = , 7 , F 5 = COSn, F f ~ 



/ ; — ~ > ™" "»J ~u / ; — r 5 

V 1 + sin u v 1 + sin u 



we have 



o 9 n n / / „ — cos u(2 + sin 2 u) \ 

ds^ = dx^ + (V + c?2 2 = V 1 + sin 2 ut> H ^ ^ L w du z + dw 2 + dw z 

V 1 + sin z n / 

Example 8. The negative constant curvature metric 

2 dx 2 + dy 2 + dz 2 (u + v) 2 du 2 + dv 2 + dw 2 



ds 

for 



'W 2 



x = (u + v ) sin it + cos it, y = —(u + v) cosn + sinn, z = w, 

leads to nonregular separation for Hamilton-Jacobi equations, and in the quantum case, to 
nonregular i?-separation for velocity dependent potentials [7]. 



9 Conclusions 

We have demonstrated that the characterization of symmetry related solutions modulo a side 
condition for Hamilton- Jacobi and Helmholtz or Schrodinger equations coincides with maximal 
nonregular separation of variables for which there is a generalized Stackel matrix with one 
arbitrary column. We have also shown that these systems can be characterized geometrically, 
i.e., in a coordinate- free manner. We have demonstrated that there is a structure theory for this 
type of separation; it is not just a collection of examples. This allows us to obtain new separable 
solutions for these equations that cannot be obtained by the standard (regular) methods of 
separation of variables. 

Our work leads to additional questions and possible extensions. 

• Eisenhart showed that the Robertson condition for Helmholtz separabilty of a Hamilton- 
Jacobi regular separable coordinate system was the vanishing of the off-diagonal elements 
of the Ricci tensor in these coordinates [12, 14]. Is there a corresponding geometrical 
interpretation of the obstruction problem for nonregular separation? 

• Is there a physical interpretation of the side conditions? 

• For regular separation of Hamilton-Jacobi and Helmholtz equations all separable systems 
on iV-dimensional Euclidean spaces and A^-spheres are known [15]. Can a similar clas- 
sification of separation with a side condition be carried out? 

• Much remains to be done in finding explicit, physically interesting, new solutions via this 
method. 

• It is clear that separation with a side condition is also appropriate for heat, Laplace and 
wave equations, where the new solutions are likely to be more interesting physically. This 
should be done and a structure theory worked out. 
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